arXiv:2205.03726v1 [astro-ph.SR] 7 May 2022 


Possible Solution to the Triple Alpha 
Fine-Tuning Problem: Spallation Reactions 
during Planet Formation 


Fred C Adams* 


*Leinweber Center for Theoretical Physics, Physics Department, University of Michigan, 
Ann Arbor, MI 48109 
è Astronomy Department, University of Michigan, Ann Arbor, MI 48109 


Abstract 


Carbon is produced during the helium burning phase of sufficiently massive 
stars through the triple alpha process. The 0+ energy level of the carbon 
nucleus allows for resonant nuclear reactions, which act to greatly increase 
the carbon yields compared to the non-resonant case. Many authors have 
argued that small changes to the energy level of this resonance would lead to 
a significantly lower carbon abundance in the universe, and this sensitivity 
is often considered an example of fine-tuning. By considering spallation re- 
actions occuring during the process of planet formation, this paper presents 
a partial solution to this triple alpha fine-tuning problem. Young stellar ob- 
jects generate substantial luminosities of particle radiation (cosmic rays) that 
can drive nuclear reactions through spallation. If the standard triple alpha 
process is inoperative, stars tend to synthesize oxygen (and other alpha ele- 
ments) rather than carbon. Cosmic rays can interact with oxygen nuclei to 
produce carbon while planets are forming. The resulting carbon abundances 
are significant, but much smaller than those observed in our universe. How- 
ever, for a range of conditions — as delineated herein — spallation reactions 
can result in carbon-to-oxygen ratios roughly comparable to those found on 
Earth and thereby obviate the triple alpha fine-tuning problem. 
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1. Introduction 


The laws of physics, as realized in our universe, are characterized by a 
collection of fundamental constants that set the strength of the forces and 
the masses of particles; additional cosmological quantities specify the inven- 
tory of the universe [I] 2]. All of these parameters have the proper values to 
allow for a long-lived universe that can develop a range of cosmic structures, 
including galaxies, stars, planets, complex nuclei, and even life. A great 
deal of previous work has considered the possibility that the fundamental 
constants could have different values, and has explored the degree to which 
they can vary and still allow the universe to be habitable [3] [4] {5} {6} [7] {8} 19]. 
Within this enterprise, an important issue is the production of carbon, an 
essential ingredient for the development of life, at least in terrestrial forms. 
Variations in the fundamental constants can lead to corresponding variations 
in the resonance levels of the carbon nucleus, and can significantly change 
the carbon production rate during stellar nucleosynthesis. The resulting car- 
bon abundance depends rather sensitively on the carbon resonance structure 
and hence on the underlying fundamental constants. This sensitivity is often 
cited as an example of fine-tuning, where small changes to the fundamental 
constants could lead to large changes in carbon production and lead to a 
lifeless universe. This paper explores an alternate channel of carbon produc- 
tion through spallation reactions that take place during the planet formation 
process and thus offers a way to avoid the triple alpha fine-tuning problem. 

Carbon is synthesized in stars during their helium burning phases and oc- 
curs through the triple alpha process [10] [I]. The expected reaction, where 
two *He nuclei are forged into a 8Be nucleus, is compromised because the ®Be 
nucleus is unstable, with a half-life of only ty/2 ~ 10716 sec. Even with this 
short half-life, however, the stellar core builds up a small transient population 
of ŝBe, which can then interact with additional alpha particles to produce 
carbon [12]. The latter reaction occurs rapidly enough to account for the 
carbon inventory of our univese, but requires a resonance in the carbon nu- 
cleus , where the resonance energy is 7.644 MeV [14], and corresponds to 
a Ot nuclear state of the carbon nucleus. If the energy level of the resonance 
were higher (lower), then the reaction rate could be slower (faster), resulting 
in different carbon yields. 

Previous work has explored the sensitivity of carbon production to the 
specific value of the carbon resonance energy level. These studies show that 
the energy level of the resonance is sensitive to the values of the fundamental 


constants [15] [6] and that the carbon yields in stars are sensitive to the 
value of the resonance [19]. If the resonance energy is higher, with 
(AER) > 0, then carbon yields from massive stars decrease. In contrast, 
lower resonance levels with (AER) < 0 lead to increased carbon production. 
The range of resonance energy that allows for net carbon production is ap- 
proximately -300 keV < (AER) < 500 keV, or a total span of ~ 800 keV 
[19]. This range should be compared to the spacing of energy levels in the 
carbon nucleus, about 3 MeV, which implies roughly a one in four chance of 
the resonance energy falling in the proper range to allow for carbon produc- 
tion. As another point of comparison, the ŝBe nucleus only fails to be stable 
by 92 keV. As a result, the changes to the energy levels of nuclei required to 
compromise the triple alpha process are almost an order of magnitude larger 
than the changes needed to make Be stable. The latter possibility allows 
for carbon production without any need for the triple alpha process [20]. 

The above discussion applies to the net amount of carbon produced in 
stars. For the parameter choices that lead to low carbon yields, the stars 
still process alpha particles into carbon. The issue is that the nuclear burn- 
ing temperature required for carbon production becomes large enough that 
newly synthesized carbon is promptly processed into oxygen (e.g., through 
the reaction '*C + *He — 160), so that relatively little carbon remains. For 
even more extreme parameter choices, much of the oxygen can be processed 
into larger alpha elements such as neon, magnesium, and silicon. As a result, 
when stars fail to produce substantial amounts of carbon, they leave behind 
oxygen and other alpha elements instead. These nuclei can be broken down 
to produce carbon through spallation reactions. 

Note that this paper implicitly considers variations in the energy levels 
of carbon of order 100 to 500 keV. Although the binding energies of carbon, 
oxygen, and other alpha elements are expected to also change by roughly 
comparable increments, such changes are small in relative terms. Recall that 
the binding energy of carbon is about 92 MeV, and that for oxygen is about 
124 MeV [I0]. As a result, changes to the binding energies are small, even 
though large changes to the carbon production rates can be realized. 

By considering spallation reactions that convert oxygen into carbon dur- 
ing the process of planet formation, this paper constructs an alternate mech- 
anism to circumvent the triple alpha fine-tuning problem. Planets form in 
the circumstellar disks associated with young stellar objects during the first 
~ 1 — 10 Myr of evolution. During this time span, the stellar hosts pro- 
duce copious amounts of particle radiation — cosmic rays — which have the 


potential to change the nuclear composition of the planet-forming material 
through spallation. In the model explored here, strong magnetic fields cou- 
ple to the disk at radii r ~ ry ~ 0.1 AU, where reconnection events lead to 
particle acceleration and large fluxes of cosmic rays (Section [2). This parti- 
cle radiation can drive spallation reactions that convert some fraction of the 
oxygen nuclei into carbon (Section B). Because the cosmic rays originate near 
the inner edge of the disk, planets forming with shorter orbital periods tend 
to have higher carbon abundances. For small stars, planets in the habitable 
zone can have carbon to oxygen ratios comparable to that of Earth (Section 
[4). As summarized in Section |5| this scenario allows for the development of 
potentially habitable planets, even in universes where the triple alpha process 
cannot produce high enough carbon yields through stellar nucleosynthesis. 


2. Production of Particle Radiation during Early Stellar Evolution 


During their early formative phases, young stars are accompanied by 
circumstellar disks that provide the initial conditions for the planet formation 
process. These disks are truncated near the host star due to the strong 
magnetic fields generated within the star. The magnetic truncation radius 
rx is determined by the balance between the inward pressure due to mass 
accretion through the disk and the outward pressure due to the fields, and 
can be written in the form 


BR y” 
a k a 


where M is the mass accretion rate and B, is the field strength on the stellar 
surface. Other parameters include the stellar radius R, and the stellar mass 
M.. Finally, 7 is a dimensionless parameter that is expected to be of order 
unity, and its value varies with the model details [22]. For typical pre- 
main-sequence stars, surface field strengths fall in the range B, ~ 1 — 2 
kilogauss, so that truncation radii ry ~ 0.1 AU (note that all of the physical 
quantities appearing in equation |1| have been measured for a large number 
of pre-main-sequence stars ). 

In the astrophysical model under consideration [24| [25], particle acceler- 
ation occurs due to magnetic reconnection events that take place near the 
truncation radius of equation (ip. The strong magnetic field lines tend to 
make the star co-rotate with the inner edge of the disk (at rx), but fluc- 
tuations in the mass accretion rate, the magnetic field strength, and other 
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parameters require the field lines to continually adjust through reconnection 
events. This magnetic activity, in turn, leads to particle acceleration. This 
process is essentially a scaled-up version of the mechanism that produces 
cosmic radiation in the Sun [26]. The luminosity L, in energetic particles is 
expected to be comparable to, but somewhat less than, the luminosity Ly 
emitted in X-rays [27]. The latter luminosity is observed to be a substantial 
fraction of the photon luminosity L,, with both temporal and source-to- 
source variations [29]. As a result, the power contributions from young 
stellar objects are ordered according to the relation 


Ly ~ 107*Lx ~ 107*E,. (2) 


Note that pre-main-sequence stars are brighter than their main sequence 
counterparts, and that low mass stars evolve more slowly than larger stars 
during PMS contraction. As a result, the typical stellar luminosity is given 
by L, ~ 1Lo during the epoch of planet formation (stellar ages t < 10 Myr). 

Cosmic rays generally have a power-law energy distribution |30} [BI], which 
can be written in the form 


eae) =9 (2), (3) 


where the distribution is normalized over the range Ey < E < oo. For cosmic 
radiation generated in the Sun through gradual flares, the power-law index 
p ~ 2.7 [B2]. For particle radiation generated through implusive flares, the 
distribution is steeper with p ~ 3.5. For the sake of definiteness, we use p = 3 
for the estimates of this paper (for completeness, note that the distribution 
for galactic cosmic rays has energy dependence ~ E~?" at high energies, 
which is roughly similar). The luminosity relations of equation (2) correspond 
to cosmic rays with energies E > Ey = 10 MeV. As outlined below, the 
energy threshold E;n for the nuclear reactions of interest is comparable to, 
but somewhat larger than, the scale Ep. For a given energy spectrum, the 
number of cosmic rays produced per unit time Nor with energies E > Ein 
is related to the luminosity L, according to 


Niog= PZB Ia ( Ba)" _ by (Ba) w 
p— 1 Eo Vy 2Eo \ En] ’ 


where the final equality assumes the index p = 3. Here, L, is the total 
luminosity in cosmic rays over the energy range E > Eo, so that the factor 


(Eo/E;,)? corrects for the inclusion of only those particles with energy E > 
Ly, > Eo. 

In this scenario, the cosmic radiation is generated within an annulus in 
the disk roughly centered on the truncation radius ry from equation (1). 
As a result, the number flux (fluence) of cosmic rays that is generated with 
energy E > Ep is given by 


N, L,E ae 
Dor & = PO 2x 10° (=) cm~? s7}, (5) 


= 2 p2 
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where the numerical estimate uses Ey = 10 MeV, L, = 10-*Lo, and ry = 0.1 
AU. 

Equation specifies the number flux in the vicinity of the truncation ra- 
dius r ~ ry. The subsequent propagation of the particles from their point of 
origin is complicated. Since the magnetic field lines are continually wrapped 
up by the shear flow from the disk, and then reconnected, the corresponding 
field structure is expected to be complex and chaotic in the region where the 
cosmic rays are accelerated. 

In this setting, the disk scale height h provides an important length scale 
in the problem, where h ~ rx /20 (see, e.g., [24]). The cosmic rays propagate 
by moving along the field lines, with a typical length scale h, and perpendicu- 
lar to the field lines, with typical length scale set by the magnetic gyroradius 
rg = ymcv/qB. For typical energy (10 MeV — 10 GeV) and field strength 
(B ~ 1 gauss), rg < h, so that the particles are tied to the magnetic field 
and travel along the highly convoluted field lines. In the absence of attenu- 
ation, cosmic ray propagation can thus be modeled as a random walk with 
step length h (note that attenuation is considered in the following section). 
In order for cosmic rays to leave the annulus where they are generated, they 
must travel a total distance d ~ ry = VN h, where N is the required number 
of steps. As a result, the cosmic rays must propagate across the disk scale 
height many times before escaping, i.e., M = (rx/h)? ~ 400. The total 
distance traveled is given by Mh = r&/h ~ 20rx. With such a long path 
length, a large fraction of the cosmic rays will be absorbed within the annulus 
where they are generated. 


3. Spallation of Oxygen to make Carbon in Circumstellar Disks 


The current working scenario for rocky planet formation can be summa- 
rized as follows. The raw materials for planet formation begin in the form of 


dust grains. In the interstellar medium, a substantial fraction of the metals 
(elements beyond helium) are locked up in these entities. Although they have 
a distribution of sizes, the grains begin with a typical size bọ ~ 0.1 — lum. 
Within the circumstellar disks that form planets, the dust grains grow into 
larger bodies, sometimes known as pebbles, with sizes b, ~ 0.1 — 1 cm. The 
pebbles sink to the midplane of the disk, experience a streaming instability, 
and thereby produce rocky planetesimals with sizes bg ~ 1 — 100 km 
(alternately, planetesimals can be produced through a gravitational instabil- 
ity [34]). These planetesimals subsequently accumulate into rocky terrestrial 
planets. 

In our galaxy the original dust grains, and hence the rocks they grow 
into, are generally composed of both silicates and graphite, with admixtures 
of many other minerals . For the case under consideration here, with little 
carbon produced in stars, the grains are expected to be primarily silicates and 
thus contain a large mass fraction of oxygen. For example, typical compounds 
found in solar system rocks include SiO and Fe2O3 (see [36]). 

If the rocky entities are exposed to a flux of cosmic radiation, then spal- 
lation can produce carbon through a number of possible reactions, e.g., 


6O+p—> PC+3p+2n and “O+p— °C+pta, (6) 


along with many others. In this case, oxygen-16 represents the target nucleus. 
For common reaction products — such as carbon — the cross sections for such 
spallation reactions are generally of order og = 100 mb (40) [41]. 
For spallation reactions that convert oxygen to carbon, the data show that 
there is an effective threshold energy at E = Ey, ~ 15 MeV, where the cross 
section rises steeply as a function of energy, and then varies slowly for larger 
energies. For the estimates of this paper, we adopt the simple approximation 
where o(E) = oo = 100 mb = 107% cm? for E > E;n and o(E) = 0 otherwise. 

Carbon production is maximized if the target oxygen nuclei are exposed 
to the full flux of cosmic rays, i.e., in the optically thin limit. In practice, this 
condition holds if the rocky bodies are in the form of pebbles (b, ~ 1 cm), 
whereas planetesimals (be ~ 1 — 100 km) are optically thick. On the other 
hand, if the bodies are too small (bj ~ 1 um), they remain coupled to the 
gas and are swept into the star along with the gas in the accretion flow. We 
thus consider intermediate sizes. These centimeter-sized bodies are expected 
to be decoupled from the gas inside the truncation radius r < rx, where the 
column density of the gas is low [24]. 


With the above specifications, the rate [ at which a given target oxygen 
nucleus interacts to make carbon is given by 


L OoLo 
T= Porla) = oar? E . 


(7) 


In the limit where the target nuclei are exposed to an unattenuated cosmic 
ray flux for an exposure time (At), the time evolution of the carbon/oxygen 
inventory is given by 


4 Xc 
3 Xo16 


L,ooEo(At) 


2 p2 
2rrg Ep, 


[C/O] = 


= exp[I'(At)] — 1 #T(At) = ~ 0.0025, (8) 


where [C/O] is the ratio of the number of carbon to oxygen nuclei and Xx 
is the corresponding mass fraction of the nuclear species (note that we also 
assume that carbon is the most common product of the spallation reactions). 
The final approximate equality holds in the expected limit where only a small 
fraction of the oxygen is converted. The numerical estimate assumes typical 
values for the parameters, i.e., Lp = 10~-*Lo, co = 100 mb, Eo = 10 MeV, 
Em = 15 MeV, rx = 0.1 AU, and (At) = 1 Myr. 

Note that the above estimate only includes reactions for carbon produc- 
tion where oxygen nuclei are the target, but other targets are possible. In 
the scenario under consideration here, we are assuming that carbon is rare 
because most of the carbon produced in stars is immediately synthesized into 
oxygen. The expected oxygen abundance is thus given approximately by the 
combined mass fraction of carbon and oxygen in our solar system, or about 
~ 2/3 of the total mass fraction of metals [86]. The remaining ~ 1/3 of 
the metal mass fraction provide additional targets for carbon production (as 
well as additional attenuation — see below). Since the corresponding reaction 
cross sections are lower, however, the correction for non-oxygen targets is 
expected to be modest. In any case, the neglect of alternate targets renders 
the estimates of this paper as lower limits. 

The expected exposure time (At) is subject to several constraints. Cos- 
mic ray production occurs due to reconnection events taking place near the 
truncation radii of circumstellar disks, so the disk lifetime limits the exposure 
time (at least for the most energetic phases). Observed circumstellar disks 
have a distribution of lifetimes that can be modeled as an exponential 
with half-life ~ 3 Myr. As a result, only ~ 10% of disks live as long as 10 Myr, 
which thus provides an approximate upper limit on (At). As outlined below, 
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another constraint arises from the growth of the rocky material containing 
the oxygen targets. While the rocks are small, they are exposed to the full 
flux of cosmic rays, but shielding occurs once the rocks have grown larger 
than b ~ 1 — 10 cm. The time at which rocky bodies cross this size thresh- 
old, typically a few Myr [43], provides another constraint on the effective 
exposure time. Calculations that explain the observed/inferred abundances 
of short lived radioactive nuclei in our solar system lead to similar time scale 
estimates [44]. In order for spallation to simultaneously produce the correct 
isotopic ratios for 7°Al/?7Al, “4Ca/*Ca, *23Mn/*Mn, and °*Nb/*?Nb, the 
irradiation time must be of order ~ 1 Myr (see also [24] |25]). 

Since the target nuclei of interest tend to be locked up in rocky bodies, 
maximal exposure occurs only for sufficiently small entities. When cosmic 
rays enter rock, they lose energy primarily through electromagnetic interac- 
tions with electrons, rather than through nuclear reactions. For example, 
cosmic rays with energy E = 10 MeV lose their energy after passing through 
a column density of only ©, ~ 0.3 g cm~? [45] 46], with larger stopping 
columns for particles of higher energy. In addition to energy loss, the prop- 
agation of cosmic rays through rock is complicated by the emission of sec- 
ondary particles, such as neutrons and electrons, which can propagate longer 
distances (although their impact is small in the present application). As a 
result, one must consider the full energy distribution of cosmic rays and per- 
form detailed simulations of the process. For example, for the case of solar 
cosmic rays impinging on the Martian surface, the isotope production rate 
for carbon decreases with depth z with an nearly exponential fall-off [46]. 
Motivated by this finding we write the production rate in the approximate 
form 

y = wexpl-2/C], (9) 


where we expect the effective scale height ¢ ~ 1 cm {45} [46] (and where 
this approximation is valid for solar cosmic rays). The corresponding rate 
coefficient yo is calculated for an incident flux of solar cosmic rays in the 
energy range E > 1 MeV. If we rescale the coefficient to include only cosmic 
rays with energy E > Ey = 10 MeV, and define a benchmark number flux of 
protons ®; = 1 cm~? s~t, then we find the value yo ~ 0.003(®cR/®1) atom 
cm? s71, where ®gp is the number flux of cosmic rays with energy E > 10 
MeV. Using these results, let us consider a cosmic ray flux cpr incident on 
a Slab of rocky material with thickness b. The rate at which oxygen nuclei 


are converted into carbon is given by the expression 


16m 2) a 
Tr = (|? } = [1 —exp(—6 ; 10 
p= (FER) $i — exp(-/0) (10) 
where pr = 2.5 g cm~® is the density of the rock. In the limit of small 
thickness b < ¢, the function on the right approaches unity, and the rate 
becomes 


ald <6) = (10 st) (S28) (11) 
1 

This rate p(b < Ç) can be compared to the optically thin rate I = (a) ®cp 
from equation (7). The rates are comparable in this limit, as expected. For 
completeness, note that the optically thick rate from equation in the 
optically thin limit is equal to the optically thin rate from equation 
(7) for a specific choice of oxygen abundance (i.e., the rate coefficient yo 
depends on Xo16). In any case, rocky bodies with size b < Ç can be treated 
using the optically thin expression. As a working approximation, we assume 
that the net effect of rocky material agregating into larger bodies of size b 
is to reduce the effective rate of converting oxygen into carbon by the factor 
~ (¢/b)(1 — exp[b/¢]). As a result, the largest rocks that remain optically 
thin have sizes less than about b ~ C ~ 1 cm. 

For larger rocks, b > ¢, only the outer layer is exposed to the cosmic 
ray flux. Within the region of disk where the cosmic rays are generated, 
strong but tangled magnetic fields confine the particle radiation, which must 
propagate diffusively. As a result, the cosmic radiation is expected to be 
nearly isotropic. For larger rocks, assuming spherical geometry, the fraction 
F of their volume that will experience spallation reactions is given by 


ra($)-3(S\ +(%) =, (12) 


where the final expression corresponds to the limit b > ¢. 

The discussion thus far assumes that the target oxygen nuclei, or rocky 
bodies that contain them, are exposed to the full flux of cosmic rays. If the 
total amount of solid material in the reconnection region is large enough, 
then the cosmic ray flux can be attenuated by background before reaching 
the rocks of interest. We can address this issue by estimating the optical 
depth of the background field of rocky material. Let the total mass M, in 
solids be spread out over the reconnection region, which is an annulus with 


10 


radius comparable to the truncation radius rx. The column density of solid 


material is thus given by 
Ms 


D, es. 13 
= (13) 


The solid material is primarily contained in rocky bodies, which can be 
considered as partially optically thick. We further consider the rocks to have 
a given radius b. As outlined above, the nuclei within the rocks are exposed 
to only a fraction F(b) of the total cosmic ray flux, whereas the correspoding 
fraction [1 — F(b)] is lost (primarily through electromagnetic interactions) | 
The optical depth of the reconnection region can then be written in the form 


M, 3 200 


[1 — F(b)] = (elem) |! =F (0) |), (14) 


Toack X mre, 4pb 
where the numerical estimate uses benchmark values M, = 1Mẹ and ry = 
0.1 AU. In the limit of large bodies b >> ¢, the fraction F — 0, so that all 
of the cosmic rays striking a given rocky body are lost. In this case, the 
optical depth Tack is determined by the surface (number) density of rocky 
bodies times their geometrical cross section. In the opposite limit of small 
rocks with b < ¢, the fraction F > 1 but 1— F = O(b/C) so that the optical 
depth becomes independent of the rock size b (e.g., see |48]}). 

The optical depth Toack and the factor F jointly determine the fraction 
fer of the total cosmic ray flux seen by the typical target nucleus. If we 
assume a simple one-dimensional geometry for cosmic ray propagation across 
the reconnection region, then the nuclei within the rocks are exposed to a 
fraction of the original cosmic ray flux given by 

1 _ 
fer = [1 =e Pach | F(b) ` (15) 


Tback 


The first factor takes into account the loss of cosmic rays due to the back- 
ground sea of rocky bodies and the second factor represents the fraction of 
nuclei within a given rock that are exposed to cosmic rays. 

The survival fraction fer depends on both the total amount of mass Ms 
in the reconnection region and the size b of the individual rocks. For solid 
masses M, < 1Mọẹ, the optical depth Ta. < 1 and the survival fraction 


'Note that the particles are still present, but they are slowed down below the energy 
threshold so that they are no longer useful for spallation reactions. 
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fe = F. If, in addition, the rocks are small so that b<1 cm, both fer and 
F are of order unity so that the nuclei are exposed to nearly the full flux 
of cosmic rays. In this limit, the carbon-to-oxygen ratio evolves with time 
according to equation (8). In general, however, the interaction rate I’ (see 
also equation m) must be reduced by the factor fer. For example, in the 
case of large total mass M, Z 1Mg, the optical depth is expected to be large 
(of order ~ 100), so that equation reduces to fer ~ F(b)/Tback- Since 
the rocks are (most likely) in the form of pebbles with b ~ 1 cm, the fraction 
F ~ 1, so that the reduction factor fer ~ 1/Tback- 

For completeness, we note that cosmic rays can also lose energy through 
interactions with gas in the reconnection region. However, such losses are 
already incorporated into this treatment, since the cosmic ray luminosity 
Lp, as defined here, corresponds to those particles that are emitted, i.e., the 
particles that survive propagation over distances comparable to the inner disk 
size. Moreover, the column density of gas in the region r < rx is expected to 
be small, ~ 1074 g cm~?, based on soft X-ray observations of young stellar 
objects [24]. For the estimates of this paper, we take into account possible 
attenutation by considering smaller values of the original (net) cosmic ray 
luminosity Lp. 

Figure |1| illustrates the possible time evolution of the carbon to oxygen 
ratio [C/O] for various levels of attenuation, corresponding to various total 
mass in solids in the reconnection region. If the rocky material resides in 
pebbles (b ~ 1 cm) and experiences the full unattenuated flux of cosmic 
rays, then the [C/O] ratio evolves with exposure time according to equation 
(8). as shown by the solid black curve. The remaining curves show the time 
evolution of the [C/O] ratio for cases where the cosmic ray flux is attenuated, 
with the total optical depth given by log;,7 = 0.5 (blue), 1.0 (cyan), 1.5 
(magenta), and 2.0 (green). These optical depth values correspond to total 
masses in solid material of M,/Mẹ = 1073/2, 10-1, 1071/2, and 1.0. The 
rocky material, which initially resides in small pebbles, eventually transitions 
into planetesimals with sizes R = 10-100 km. At this point, the the rate of 
transforming oxygen into carbon decreases abruptly by a factor of ~ 104—105. 
As a result, the [C/O] ratio effectively freezes out at the transition time, i.e., 
the transition time defines the effective exposure time. 

We can summarize this scenario for carbon production as follows: The 
rocky material is expected to reside in the form of pebbles with size b ~ 1 cm 
and the total mass in the reconnection region is comparable to that of a small 
rocky planet (so that M, ~ 0.1 — 1Mẹ). Finally, the exposure time falls in 
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Figure 1: Time evolution of the carbon to oxygen ratio [C/O] for varying levels of at- 
tenuation. The upper black curve shows the maximum growth of [C/O] as a function of 
time for optically thin conditions and the full expected flux of cosmic rays. The remaining 
curves show the [C/O] ratio for total optical depths given by Tback ~ 3.16 (blue), 10.0 
(cyan), 31.5 (magenta), and 100 (green). For the expected case where the rocky material 
resides in pebbles (b ~ 1 cm), these optical depths correspond to total mass in solids 
M,/Ma 7% 0.0316, 0.1, 0.316, and 1.0, respectively. 
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the range (At) = 1—10 Myr, consistent with observed disk lifetimes. Under 
these conditions, the optical depth of the background is large, Tpack >> 1, so 
that the survival fraction fer oc Til, « My. Moreover, only a small fraction 
of the total number of oxygen nuclei will be processed into carbon. In this 
regime of parameter space, the carbon to oxygen ratio can be written in the 


rn teyoy oan (88) (27. r 


4. Planet Formation and Habitable Zones 


The previous sections show that cosmic rays are accelerated near the 
truncation radius ry, are confined to nearby regions, and can interact with 
oxygen nuclei to make carbon. This section briefly considers the conditions 
necessary for the resulting carbon to be incorporated into a potentially hab- 
itable planet. 

Even under the most favorable conditions, the carbon abundance is quite 
small, only a fraction of the total mass of the planet (which is made from the 
rocky raw material enriched via spallation). However, the Earth itself also 
has relatively little carbon, in contrast to its large abundance in both the 
Sun and the Galaxy. Even the carbon abundance for Earth’s crust has been 
difficult to measure, with estimates of the carbon to oxygen ratio falling in 
the range [C/O] ~ 0.0002 — 0.01 [49] 50, 51) (52) (53) [54]. Estimates for the bulk 
composition of Earth are more uncertain, but fall in a similar range. Finally, 
we note that oxygen makes up ~ 30% of the Earth’s mass, so the mass 
fraction of carbon is smaller than the carbon to oxygen ratio, say Xcj2 ~ 
0.0001 — 0.003. Significantly, the [C/O] ratios that are produced through 
spallation reactions (Figure can be comparable to (and under extreme 
conditions larger than) the values estimated for Earth. However, the largest 
[C/O] ratios only arise for relatively small total mass M, in solids. To form 
planets with M, = M, = 1Mz and an exposure time of 10 Myr, we find 
[C/O] ~ 0.0003, which is near the low end of the range inferred for Earth. 

The most straightforward scenario is for the planet to form at the annulus 
centered on rx where the cosmic rays are generated. In this case, the carbon 
to oxygen ratio in the rocky material will sample the values illustrated in 
Figure |1| i.e., [C/O] ~ 1075 — 10°. From this birth scenario, the resulting 
planetary orbit has a semimajor axis near a ~ ry ~ 0.1 AU. In order for the 
planet to receive the same stellar flux as Earth, the star must be less luminous 
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than the sun by a factor of 100, which corresponds to a main sequence star 
with mass M, = 0.25Mo. 

Figure |2} shows that rocky planets readily form near the location of the 
magnetic truncation radius where cosmic rays are generated. The blue his- 
togram in the figure depicts the semimajor axes for a collection of exoplanets 
detected by the Kepler mission [55]. This sample includes the planets with 
inferred masses Mp < TMọẹ that are found in systems with three or more 
planets. Planets with masses below this threshold exhibit a primarily rocky 
composition, whereas larger planets tend to have extended atmospheres (al- 
though substantial variation exists). Only planets in multi-planet systems 
were used, as those tend to have nearly circular orbits and make up a uni- 
form sample|?| The green vertical lines show the locations of the magnetic 
truncation radii rx for a collection of pre-main-sequence stars with measured 
magnetic field strengths [23]. The annulus where cosmic rays are accelerated 
through reconnection and flaring activity is expected to extend mostly inward 
from rx [25]; to illustrate the extent of the enrichment region, the vertical 
cyan lines show locations corresponding to rx /2 for the same sample. 

Note that the present-day locations of the planets (as shown in the figure) 
are not necessarily where they were formed. If rocky planets form at larger 
radial distances from their host stars, and subsequently migrate inward, then 
they would not be enriched in carbon through the mechanism considered here. 

The stellar hosts for the planets shown in Figure |2| have masses in the 
range M, ~ 0.10 — 1.35Mọo, whereas the stellar masses for the rx sample 
have M, ~ 0.3 — 2 Mo. The stellar masses in both data sets thus span the 
same range and include the low mass stars of interest (where the habitable 
zones are roughly coincident with the disk truncation locations). Moreover, 
neither the planet masses M, or the truncation radii rx show any obvious 
trends with stellar mass. 

Since we are interested in potentially habitable planets around smaller 
stars, it is important to note that M stars are generally more active than 
solar-type stars and are subject to intense flaring activity. Some authors 
have speculated that such energetic activity could be detrimental for plane- 
tary habitability. Although the conditions required for life remain unknown, 
this complication should be noted (see the discussion of and references 


?Keep in mind that for many planets the mass is inferred from measurement of the 
planetary radius, i.e., the mass is not directly determined. 
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Figure 2: Histogram showing semimajor axes of observed rocky planets (thick blue curve) 
and the inferred positions of the magnetic truncation radius (vertical lines). The planet 
sample includes detected planets from the Kepler mission that have estimated masses 
Mp < 7Me@ and are found in systems containing 3 or more planets. The green lines show 
the locations of the magnetic truncation radii rx for measured properties of pre-main- 
sequence stars [23] according to equation (ip. The cyan lines show locations corresponding 
to half of the inferred values (rx/2). The region where cosmic rays are generated, as 
delineated by the collection of vertical lines, overlaps with the observed locations of rocky 
planets. 
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therein). On the other hand, the enhanced flaring activity of these small 
stars could result in enhanced spallation rates and somewhat larger carbon 
yields. 


5. Conclusion 


This paper presents a new possible solution to the triple alpha fine-tuning 
problem by considering how spallation changes the chemical composition of 
rocky planets forming in tight orbits around low-mass stars. If the funda- 
mental constants of nature had different values, so that the Hoyle resonance 
in the carbon nucleus had a different energy level, then carbon production 
could be greatly reduced. This paper shows that even in universes that fail 
to produce carbon by stellar nucleosynthesis, some carbon is produced by 
spallation outside of stars. This argument can be summarized as follows: 

If stars fail to produce substantial carbon yields, the nuclei are processed 
instead into oxygen and other alpha elements. These larger nuclei can be 
broken down into carbon by cosmic radiation through spallation reactions, 
thereby providing an alternate source of carbon. 

Cosmic ray acceleration takes place near the magnetic truncation radius 
rx in the disks associated with young stellar objects. Rocky planets readily 
form in this region, which is roughly located at a ~ rx ~ 0.1 AU (Figure 2}. 
This region corresponds to the habitable zone for stellar hosts with masses 
M., ~ 0.25Mo, which are the most common stars in our universe. As a result, 
the raw materials that make up this class of habitable planets are naturally 
irradiated by local (stellar) cosmic rays and are subject to spallation. 

The amount of carbon produced depends on the timing of events that 
take place during the process of planet formation. For the expected cosmic 
ray luminosities, if the rocky material remains optically thin for a typical 
time scale of 1 Myr, then the resulting carbon to oxygen ratio [C/O] ~ 107%, 
comparable to the value inferred for Earth. However, this value is only 
realized for relatively small total masses Ms < Mẹ. Longer exposure times 
lead to higher carbon yields, whereas larger masses M, or larger rock sizes b 
lead to greater attenuation of the cosmic ray flux and hence smaller carbon 
yields (Figure|1). For example, if the total mass in solids in the reconnection 
region is comparable to Earth, then an exposure time of 10 Myr leads to 
[C/O] ~ 0.0003, near the low end of the range of inferred values for Earth. 

The results of this work have a number of implications: Because spal- 
lation can produce carbon in universes where the triple alpha process is 
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inoperative, the fine tuning of the fundamental constants necessary for a 
viable universe is less severe than previously claimed. On the other hand, 
this spallation mechanism tends to produce only modest amounts of carbon, 
[C/O] ~ 1074 — 1078, values roughly consistent with Earth abundances, but 
much smaller than that of the Sun or the universe as a whole. In addition, 
carbon production takes place within the habitable zones of smaller stars 
with masses M, ~ 0.25 Mo. Although such hosts are the most common 
stars in our universe, the habitable zone for solar-type stars lies farther away 
from the cosmic ray source, so that direct analogs of our solar system would 
remain carbon poor. At the present time, we do not know exactly how much 
carbon is required for a planet, or a universe, to be habitable. Given the 
modest levels of carbon produced via spallation, however, it is likely that 
such universes would be somewhat less habitable than our own. 

Another lesson from this consideration of spallation is that universes — 
including our own — have a number of possible channels to reach habitabil- 
ity. For example, a large number of different astrophysical sources contribute 
to the production of carbon and other elements, including supernova explo- 
sions, collisions of compact objects, red giant winds, spallation, and big bang 
nucleosynthesis (although the latter tends to make only light nuclei). If one 
channel of carbon production is unavailable, then astrophysics (often) pro- 
vides alternate mechanisms. In the present context, spallation can produce 
carbon in other astrophysical settings, including the interstellar medium and 
on the surfaces of mature planets orbiting main sequence stars (as outlined 
in Appendix A}. 

Finally, we note that spallation from locally produced cosmic rays can 
have important implications for the chemical enrichment and early evolution 
of planets in our universe. In this context, carbon is already abundant, so 
any additional carbon produced via spallation is negligible. Spallation also 
produces short-lived radioactive nuclei, such as 7°Al, “Ca, °3Mn, !°8La, 'Be, 
and others. These radionuclides provide important sources of heating and 
ionization during planet formation, and meteoritic evidence indicates that 
our solar system was enriched in these nuclei relative to cosmic abundances. 
A great deal of past work [24] has focused on radioactive enrichment 
mechanisms for solar system meteorites, at distances of order ~few AU. Since 
many rocky planets form much closer to their host stars (e.g., Figure), they 
are expected to experience much greater levels of radioactive enrichment, 
which can influence the formation and properties of such systems [48]. 
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Appendix A. Alternate Scenarios for Carbon Production 


As outlined in the main text, intense cosmic ray production is expected 
to continue only while circumstellar disks remain intact, typically for several 
Myr. After this time, the planets are largely formed, and the magnetic 
reconnection region at the inner disk edge is no longer present. Over longer 
spans of time, however, carbon production can take place on the surfaces of 
completed planets, or in the interstellar medium. Spallation in the interstellar 
medium leads to [C/O] ~ 1077 over a Hubble time [9]. This estimate assumes 
present-day cosmic ray fluxes, although they could be somewhat larger in the 
past (since cosmic rays are produced through supernovae, which trace the star 
formation rate, which was larger in the past). The flux could also be larger in 
other galaxies in other universes. In any case, the interstellar medium levels 
of [C/O] are much smaller than the abundances realized in the reconnection 
regions of circumstellar disks. Continued irradiation of mature planets leads 
to similarly low [C/O] levels, as considered in this section. 

After disk dissipation, when planets are fully formed, the host stars con- 
tinue to generate cosmic radiation near their stellar surfaces, with a reduced 
particle luminosity. In this context, however, we are interested in small stars, 
M dwarfs with M, ~ 0.25Mo, which are observed to have more energetic flar- 
ing activity than the Sun. As a result, the expected cosmic ray fluxes from 
M dwarfs will be significantly larger than solar values during the main se- 
quence phase (but smaller than the fluxes generated at rx during 
the pre-main-sequence phase). More specifically, estimates for the cosmic ray 
flux from M dwarfs indicate that the cosmic ray flux 6 ~ 10° cm~? s7}, for 
particles with energy E > Ey = 10 MeV and for a distance r = 0.1 AU from 
the star [59]. This estimate assumes that the M dwarf continually produces 
energetic flares (which lead to cosmic ray acceleration) and thus corresponds 
to the upper end of the possible range. With this flux, the coefficient for the 


carbon production rate in equation (9) becomes yọ ~ 300 atoms cm~? s~1t. 
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Now consider a planet orbiting a small star at 0.1 AU and assume that it 
is exposed to cosmic rays over time (At). The planet is likely to be spinning, 
but at any given time an area A ~% TR? of the planetary surface will be 
exposed to the flux of energetic particles. The total number of carbon nuclei 
produced through spallation takes the approximate form 


Norz © (1 R2)Cyo(At) ~ 10° ( Re) (4) (A.1) 


For comparison, the biosphere of Earth contains 400 — 800 billion tons of 
carbon, which is equivalent to about 2 — 4 x10*° carbon atoms. As a result, 
the number of carbon nuclei that can be produced by spallation over a Hub- 
ble time corresponds to less than one percent of the carbon content of our 
biosphere. With the assumed cosmic ray flux, the star emits a total number 
of cosmic rays Nor ~ 10° and the planet intercepts N, ~ 4 x 10*° over the 
time (At), so that Nor > Np > Noz. 

The earlier phase of stellar evolution can produce carbon to oxygen ratios 
up to [C/O] ~ 0.00025 for M, = 1Mg (see Figure [I). Assuming an initial 
oxygen mass fraction Xo16 = 0.3, the total carbon content for an Earth-like 
planet corresponds to N ~ 2 x 106 nuclei. This value far exceeds that total 
amount of carbon produced through the long term channel of spallation. On 
the other hand, long term exposure produces carbon at the planetary surface, 
where it could be more biologically useful. 

Looking at these results another way, the number of carbon nuclei pro- 
duced through long-term enrichment corresponds to the number produced 
through early enrichment in the upper ~ 1 cm of the planet (in the absence 
of a planetary magnetic field strong enough to shield the planetary surface 
from cosmic rays). This thickness is comparable to the scale height for spalla- 
tion (assuming Ç ~ 1 cm). As a result, long-term irradiation roughly doubles 
the carbon content of the upper layer of the planet (1 cm thick). 
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